Quantum criticality of dipolar spin chains 
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We show that a chain of Heisenberg spins interacting with long-range dipolar forces in a magnetic 
held h perpendicular to the chain exhibits a quantum critical point belonging to the two-dimensional 
Ising universality class. Within linear spin-wave theory the magnon dispersion for small momenta k 
is [A 2 + Vfrk 2 ] 1 / 2 , where A 2 oc \h— h c \ and v\ tx | lnfc|. For fields close to h c linear spin-wave theory 
breaks down and we investigate the system using density-matrix and functional renormalization 
group methods. The Ginzburg regime where non-Gaussian fluctuations are important is found to 
be rather narrow on the ordered side of the transition, and very broad on the disordered side. 
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I. INTRODUCTION 

The long-range nature and spatial anisotropy of the 
dipole-dipole interaction in quantum many-body systems 
can give rise to unconventional effects such as exotic or- 
dered phases and excitation spectra.^ The experimental 
study of these phases is now possible due to substantial 
progress in controlling the parameters of trapped ultra- 
cold atoms. While both bosonicPSI anc j fermionirP dipo- 
lar quantum gases have been realized, it remains a chal- 
lenge to design purely dipolar spin systems by localizing 
ultracold atoms or molecules with permanent magnetic 
or electric moments on an optical lattice. A promising 
strategy to obtain experimental realizations of dipolar 
magnets with localized spins uses trapped ions, which 
were recently employed to design spin Hamiltonians with 
controllable interactions between the spinsP^ 

Heisenberg magnets with dipole-dipole interactions in 
two and three dimensions have been inve stig ated the- 
oretically for more than half a centurypESl but one- 
dimensional dipolar spin chains have not received much 
attention. This may be due the fact that in condensed 
matter systems the exchange interaction is usually much 
larger than the dipole-dipole interaction, so that it has 
not been possible to realize experimentally purely dipolar 
spin chains before the emergence of the field of ultracold 
atoms. Recently several authorjSM^ pointed out that 
tunable spin chains with dipole-dipole interactions can 
be derived from two-component dipolar gases as effective 
models for the spin degrees of freedom. The investiga- 
tions presented in this work are motivated by the expec- 
tation that in the near future it will be possible to design 
purely dipolar spin chains using trapped atoms or ions 
at ultra-low temperatures. 

In a previous study of dipolar spin chains^ the long- 
range dipole-dipole interaction was truncated at the next- 
nearest neighbor, which misses the logarithmic correction 
to the spin-wave velocity discussed below. Moreover, the 
spin-wave calculations of Ref. [14] did not take into ac- 
count the tilted geometry of the classical ground state 



in the low-field phase, thus missing the quantum critical 
point which separates the tilted phase from the high-field 
phase where all spins align with the magnetic field. 

II. CLASSICAL GROUND STATE AND 
SPIN- WAVE EXPANSION 

We consider a chain of quantum spins Si of spin S in 
an external magnetic field h perpendicular to the chain 
which are coupled by both dipolar and exchange inter- 
actions. Choosing our coordinate system such that the 
chain lies along the rc-axis and the external magnetic field 
h = hz points in the z-direction, our Hamiltonian reads 

1 (j 2 
W = _ 2 ^ \ Xi - x-\ 3 : ' , ' S ''' S '' ~~ Si ' Sj ^ 

-^J^Si-Sj-h^S*, (1) 

ij i 

where sums are over the N sites x$ of a one-dimensional 
lattice with spacing a. The long-range dipolar interac- 
tion is characterized by an effective magnetic moment 
/i = gi-J-B, where g is the effective gyromagnetic factor 
and [Ib is the Bohr magneton. We assume that the spins 
are also coupled by nearest-neighbor ferromagnetic ex- 
change interactions, i.e., Jij = J > if \%i — Xj\ = a 
and — otherwise. Due to the competition between 
the Zeeman energy, which favors alignment of the spins 
along the z-axis, and the dipolar interaction, which fa- 
vors spin alignment along the ic-axis, the spins align with 
a finite tilt angle # relative to the field direction below 
a certain critical field h c , as shown in Fig. [T] The mag- 
netization points then in the direction of the unit vector 
rh = sin fix + cos-dz. The tilt angle § in the classical 
ground state (S —> oo) can be determined by replacing 
the spin operators Si in Eq. ([I]) by classical vectors Srh 
of length S, which yields the classical ground state energy 
per site 

Ho/N = -JS 2 - hS cost? - L> ^(3sin 2 7? - l)/6, (2) 
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FIG. 1. (Color online) Classical ground state of the dipolar 
spin chain defined in Eq. for h < h c . Big arrows represent 
the vectors Srh at the lattice sites (black dots). Small arrows 
represent the tilted basis vectors m and e\. The third basis 
vector e.2 — y points into the plane of the paper. 



where D is controlled by the dipole-dipole interaction 
N 
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Minimizing Hq with respect to the tilt angle $ we find 
cos $ — h/Do for h < Do, and d — for h > Do- In this 
work we show that at the critical magnetic field h c where 
the tilt angle i? vanishes the system undergoes a contin- 
uous quantum phase transition. Classically the critical 
field is h c — D Q , but for small S the value of h c is sub- 
stantially smaller than Dq due to quantum fluctuations. 
Since Hamiltonian (JlJ has no continuous spin symme- 
try, and the Z 2 symmetry Sf — > —Sf is spontaneously 
broken in the tilted phase, we expect that the quantum 
phase transition belongs to the universality class of the 
two-dimensional Ising model. 

To calculate the spin-wave spectrum, we expand the 
spin operators in the tilted basis {e.\,e%, m} shown in 
Fig. [TJ Introducing spherical basis vectors eP = e% + 
ipe 2 with p = ±, we write Si = sjrh + and = 
\ Sn=± S~ p e p . We then express the spin components in 



terms of canonical boson operators bi using the Holstein- 
Primakoff transformation, sf = S - b\bi, S+ = (5r)t = 
[25 — b[bi\ l / 2 bi. Retaining only quadratic terms in the 
bosons and Fourier transforming, bi = d^J^^e 
our bosonized spin Hamiltonian is approximated by 



H « n + [Mh + -f (b\bl k + b_ k b k ) 



(4) 



with 

A k = 
and B k 



Do 
3 



Dk 
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(3sin 2 i?-l) + J - Jk + hcosti (5) 



cos where J k = 2J5 cos(fca) and 
-Do cos(nka) 



■/?■' 



(6) 



The infinite series ^ 



oo cos(nka) 



represents the so-called 



Clausen function Ci3(fca) = ReLi3(e lfca ), where Lis(z) is 
the polylogarithm.^From the known series expansion of 
Li3(e A1 ) we obtain for \ka\ <C 1, 

D k /Do = 1 - \{ka) 2 d 1 \n{d 2 /\ka\) + 0((kaf), (7) 



where d\ = 1/[3C(3)] and d 2 = e 3 / 2 . Using a Bogoli- 
ubov transformation to diagonalize the Ha miltonian Q 
we obtain the magnon dispersion E k = W 'A 2 — |-Bfc| 2 , 
which is shown graphically in Fig. [2j For simplic- 
ity we set J = from now on. For \ka\ <C 1 and 
| h — Dq\ <C £>o the magnon dispersion is then approxi- 
mated by E k s» ^/A 2 + v 2 k 2 , where the square of the gap 
is A 2 = D (h-D ) for h > D , and A 2 = 2D (D -h) for 
h < Do- The squared velocity v\ exhibits a logarithmic 
divergence for small wave vectors, v\ = v 2 ln(l/|fca|) +Cq, 
with Vq = (D a) 2 di and Cg = |u 2 . The logarithmic di- 
vergence of v k is a unique signature of the l/|x| 3 -decay 
of the dipolar interaction in one dimension. Indeed, log- 
arithmic corrections to the dispersion are characteristic 
of one-dimensional systems with long-range interactions 
decaying as l/|x|^ with odd (3 (see Refs. [Pol and fTTl for 
the case of the Coulomb interaction, /3 = 1, and Ref. [TS] 
for the case of an arbitrary /3). 

Despite the unusual logarithmic correction to the dis- 
persion in Gaussian approximation, it seems at first 
glance that our spin-wave approach remains valid for all 
values of the magnetic field. This is not the case, how- 
ever, because in a narrow range of magnetic fields close 
to -Do the leading quantum correction to the magnetic 
moment to per site completely overwhelms the classical 
result m ss S. Retaining the leading 1/5 quantum cor- 
rection to to, we obtain 



2E k 



Because E k=0 vanishes for h D while A k=0 remains fi- 
nite for sufficiently small \h — Do\, the last term in Eq. ^ 
becomes arbitrarily large, which is also evident from the 
numerical evaluation of Eq. ^ shown in Fig. [3] (a). 

Obviously, for any S < oo there is a range of magnetic 
fields where the tilt angle is small and the leading quan- 
tum correction to the magnetization is larger than the 
classical result to S. In this regime our simple spin- 
wave approach breaks down and we need more sophisti- 
cated methods to investigate the behavior of the system. 
We have studied the quantum critical regime \h— h c \ < h c 
using both the numerical density-matrix renormalization 
group (DMRGJP^ and the analytical functional renormal- 
ization group (FRG^P methods. 



III. DMRG APPROACH 

Using the DMRG, we have calculated the ground state 
and its magnetization for systems of up to N — 128 sites 
and a varying bulk magnetic field h in the z-direction. 
In the simulations, we kept up to 320 density-matrix 
eigenvalues, leading to a maximum discarded weight of 
10~ 12 . It should be noted that DMRG simulations of 
the dipolar spin chains are numerically demanding for 
several reasons. First, the long-range nature of the in- 
teraction leads to increased correlations between differ- 
ent parts of the system that have to be encoded in the 
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FIG. 2. (Color online) Graph of the spin- wave dispersion Ek for J = and ft < Do (left), and for ft > Do (right). In the middle 
we show the gap A as a function of (ft — Dq)/Dq. 
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FIG. 3. (Color online) (a) Spin-wave result for the one-loop 
corrected total magnetic moment m given in Eq. Q as a 
function of ft/ft c for S = 1/2. The curves are for J = and 
h c = Dq. The inset shows m/S versus the logarithm of the 
reduced magnetic field \h — h c \/h c for ft<ft c (solid line) and 
for ft > h c (dashed line), (b) Corresponding DMRG results. 
The critical magnetic field ft™(oo)^0.507A) was determined 
by extrapolation of the critical fields ft™ (TV) of AT-site chains 
with open boundary conditions, as shown in the inset. 



variational ground state. Hence, the number of states 
that need to be kept is rather high compared to the 
usual spin chains with nearest-neighbor Heisenberg in- 
teraction. Second, the SU(2) spin symmetry is broken, 
so that S z is not a good quantum number. Therefore, 
the restriction of the basis states to a certain spin com- 
ponent S z , which is normally used to significantly in- 



crease the efficiency of the DMRG, is not possible. Fi- 
nally, for h<h c the ground state is two-fold degenerate 
due to the Z 2 symmetry Sf — > —Sf. Accordingly, it is 
crucial that we mix two states into the density matrix 
and target both states of the ground-state doublet. In 
the h/h c —>0 limit, the system is fully polarized in either 
the +x- or the — x-direction, but while the Z 2 symmetry 
is spontaneously broken in the infinite chain, the vari- 
ational DMRG ground state, which always describes a 
finite system, is an arbitrary linear combination of two 
states with opposite polarization directions. Since the 
mixing angle is not fixed, the measured magnetization in 
the x-direction is random. We have found that we can 
obtain reliable results by adding a small local magnetic 
field h x — 10 _10 /i 2 /a 3 on the end sites. This field explic- 
itly breaks the Z 2 symmetry and, consequently, leads to 
a unique ground state while generating an energy differ- 
ence of the order 10 _10 /z 2 /a 3 in the formerly degenerate 
doublet states. 

In finite systems, the correlation length £ is bounded 
by the system size L. Therefore, phase transitions are 
strictly possible only in infinite systems; the inverse size 
1 /L plays the role of an additional parameter that moves 
the system away from the critical point. To determine 
the critical parameters in the thermodynamic limit, the 
dependence of thermodynamic quantities on the system 
size can be investigated with finite-size scaling theory. In 
this context, one uses the notion of a pseudo-critical field 
h c (L) associated with indications of critical behavior in 
a finite system of size L. Note that, for a given system, 
multiple definitions of h c (L) that result in the correct 
value of h c in the limit L — > oo are possible. In gen- 
eral, one finds that the pseudo-critical field h c (L) scales 
as L~ 1 / l/ , where v is the correlation-length critical ex- 
ponent and £ oc \h— h c \~". The order parameter scales 
as m x (h e ,L)<xL~P/ v , with /3 the order-parameter criti- 
cal exponent. For the two-dimensional Ising universality 
class, the relevant critical exponents have the values v=\ 
and ^=1/8. 

In Fig. [3] (b) we show our numerical result for the total 
magnetic moment 
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FIG. 4. (Color online) DMRG results for the order pa- 
rameter m x = jf ~}2i{Si) as a function of [h — h' c ]/h' c for 
5 = 1/2 and different N. The main plot shows the scaled 
magnetization N^'"m x with exponents /3 = 1/8 and v = l of 
the two-dimensional Ising model versus N[h — h' c ]/h' c , with 
h' c = 0.516Do chosen to optimize the collapse of the data. The 
inset shows the raw data as a function of h/ Do. 



Since m y vanishes for the field direction selected here, 
the total magnetic moment is determined by m x and m z 
only. While the value of m x drops from about one half 
to zero in the vicinity of the pseudo-critical field h c (L), 
m z increases linearly for h < h c (L) and enters a satu- 
ration regime for h > h c (L). Consequently, m exhibits 
a minimum close to h c {L). By determining the position 
/i™(iV) of this minimum for systems of different num- 
bers of sites N—L/a, we obtain an approximate value 
/i™(oo)sa0.507£>o upon extrapolating to N — > oo. 

In contrast to the spin- wave results in Fig. [3] (a), the 
total moment is finite for any h and exhibits a large asym- 
metry, indicating that the Ginzburg regime where non- 
Gaussian fluctuations are important is very broad on the 
disordered side (h > h c ) of the transition. This corre- 
sponds to a regime where the order parameter m x is zero 
but the magnetization m z along the field direction is not 
yet saturated due to large correlations (SfS*) in the x- 
component of the spin. 

In order to ascertain whether the phase transition be- 
longs to the Ising universality class, we also investigate 
the behavior of the order parameter m x separately. In 
Fig. [4j we plot the scaled magnetization N^ v m x {h,N) 
with v=l and /3=l/8 versus N[h — h' c )/h' c for different 
values of N=L/a, where the critical field h' c is chosen to 
optimize the collapse of the data. The collapse of the 
data onto a single curve for a critical field h' c ~ 0.516£>o 
shows that our numerical results for the magnetization 
are consistent with the two-dimensional Ising universal- 
ity class. The value h' c estimated from the best collapse of 
the data and the value /i™(oo) obtained by extrapolating 
the minimum of m/S in Fig. [3] (b) differ by only 2.1%. 
We ascribe the imperfect collapse of the N—48 data to 
the fact that the dipolar interaction is long-range so that 
sufficiently large systems have to be treated before the 



true critical behavior manifests itself. 



IV. FRG APPROACH 

To investigate our model for general 5* and to verify 
that the RG flow in the vicinity of the critical point in- 
deed is consistent with the Ising universality class, we 
have studied the effect of spin-wave interactions using 
the functional renormalization groupP^ Since we are in- 
terested in the critical fluctuations, we may simplify 
the calculations using the hermitian field parametriza- 
tion of the spin-wave interactions^ in which we ex- 
press the Holstein-Primakoff bosons bk in terms of two 
canonically conjugate hermitian field operators, bk = 
[4>k\fhk + iHk/Vhk\/V2, where h k = A k + \B k \. Since 
the fluctuations of the canonical momentum 11^ remain 
gapped at the quantum critical point, we may integrate 
over the field life in Gaussian approximation and obtain 
an effective Euclidean action S^ff^] for the field 4>k de- 
scribing the critical fluctuations. The Gaussian propaga- 
tor of the 0-field is then G (K) = [lo 2 + £^] _1 , where 
K = (fc, iui) denotes both momentum and frequency, and 
the effective action S c s [4>] of the critical fluctuations can 
be expanded as 



1 1 00 

SM = ^ E GoHK^-k^k + ^ E 

K n=0,n- 

-K n fi^ ( d l) (ki... k n )4> Kl ■■■4>K. 
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X E S Xi 
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(10) 



with T the temperature. The interaction vertices 
lo (fci . . . k n ) can be expressed in terms of the Fourier 
transform Dk of the dipole-dipole interaction by expand- 
ing the spin operators in powers of the Holstein-Primakoff 
bosons bk and then setting bk — > 4>k\/hk/2. For vanish- 
ing external momenta, the interaction vertices have finite 
limits Fq™' 1 oc S ,1_ t j and in the symmetric phase (h > h c ) 
the odd vertices vanish. 

The magnon spectrum can be obtained from the poles 
of the true propagator G(K) = [Go\K) + E(iT)]- 1 , 
where 'E(K) is the irreducible self-energy of the effective 
field theory defined in Eq. (10). Since the higher order 



vertices involve increasing powers of 1/S for large S, one 
could calculate S(fsT) perturbatively. However, the cor- 
rections become arbitrarily large for h — > h c , which is 
not surprising because all higher order vertices of our 
1 + 1-dimensional field theory are relevant at the Gaus- 
sian fixed point, with canonical dimension +2. To reg- 
ulate the infrared singularities, we replace the inverse 
Gaussian propagator by Gq\(K) = Gq 1 (K) + R\(k), 
using the regulator function R\{k) proposed by LitimP^l 
In the limit of vanishing flow parameter A, the regulator 
vanishes, so that we recover our original model. 

We use a simple truncation of the formally exact hier- 
archy of FRG flow equations^! for scalar field theories of 



the type ( 10 1 to re-sum the perturbation series, expand- 
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FIG. 5. (Color online) FRG results for the critical field h ( c S) 
as a function of the spin S. For the relevant values of S > 1/2 
that might be investigated in a future DMRG study we find 



= 0.79£> , h ( c 3/2) 



0.83A), and h 



(2) 



0.85D . 



ing the flowing self-energy as 

E A (X) = E A (0) + Y A k 2 + {Zl 1 - 1 )uj 2 + 0(k 3 , w 3 ), (11) 

with flowing coupling constants Ea(0), Y\, and Z\. 
The scale-dependent long-wavelength magnon spectrum 
is then 



the logarithmic correction to the spin-wave velocity be- 
comes marginally irrelevant at the fixed point, namely 
v\ ~ [ln(A /A)]- r for A -> 0, with r 0.017, while 
c\ remains finite. This signals that scale invariance is 
eventually restored at the critical point, strengthening 
the confidence in the finite-size scaling analysis of the 
DMRG results. 

Besides investigating the behavior of the RG flow in 
the vicinity of the non-Gaussian fixed point, we may also 
use the FRG approach to estimate the critical field he 
for different values of the spin S. For a given magnetic 
field h < h c , we fine-tune the tilt angle $(h) so that the 
linear vertex TP vanishes for A — > and fix the value 
of the critical field h c where the magnetization m drops 
more steeply, corresponding to the limit i?(/i c ) —> 0. For 
S = 1/2 we obtain h { c 1/2) « 0.72D , which is significantly 
smaller than the classical result h c = Dq. However, the 
DMRG result h c « 0.51Do is even smaller, showing that 
our FRG truncation probably still misses some sizable 
effects of quantum fluctuations. In Fig. [5] we plot our 
FRG results for higher values of the spin. 



E A,k - A/ 



[vlHl/\ka\) 



•41 



(12) 



CONCLUSIONS 



where the renormalized squared gap is A A = Z\[A 2 + 
Sa (0)] , and the renormalized magnon velocities are given 
by v\ = Zav 2 and c A = Za[cq + Y\\- Note that 
in general the low energy expansion of the self-energy 
can also contain terms proportional to k 2 ln(l/|fca|) and 
oj 2 ln(l/|fca|), which are of the same order as the terms re- 



tained in Eq. (111. However, we find that the flowing self- 
energy Y> A (K) is analytic in K = for any finite A, due 
to the presence of the regulator function R\(k), so that 
these terms do not appear in our FRG approach. More- 
over, in our truncation we retain only the momentum- 



and frequency-independent parts T A 



(3) 



and r^ 4) of the 
three-point and four-point vertices. This truncation is 
then not sufficient to calculate the critical exponent rj for 
the anomalous dimension, which is determined by the 
frequency-dependent part of the four-point vertex. 

Analyzing the structure of the FRG flow equations in 
the vicinity of the quantum critical point, corresponding 
to a non-Gaussian fixed point of the flow equations,^ we 
find that the FRG confirms the Ising universality class 
expected from the general symmetry arguments and nu- 
merically derived within the DMRG analysis. For h = h c 
the wave function renormalization factor Z A and the 



rescaled three-point vertex cx Z'^ Z T ( ^ '/(Aa) 2 flow 
to zero in the limit A —> 0, while the rescaled four- 
point vertex 7 A 4 -* cx Z^T^/(Aa) 2 approaches a finite 
fixed-point value. In particular, the three-point vertex is 
marginally irrelevant at the fixed point, flowing asymp- 
totically as 7 A 3) - [MAo/A)]- 1 / 2 in the limit A ->• 0. 
(Here Ao ~ 1/a is the initial RG scale, corresponding 
to an ultraviolet cutoff). Importantly, we also find that 



z3/2 r (3) 



In this work, we have calculated the magnon spectrum 
and the magnetization curve of dipolar spin chains in 
a transverse magnetic field using spin-wave theory and 
renormalization group methods. We have shown that at 
a critical field h = h c the system exhibits a quantum 
critical point belonging to the two-dimensional Ising uni- 
versality class. Outside the critical regime, where the 
excitation spectrum is well described by linear spin- wave 
theory, the magnon velocity exhibits a logarithmic de- 
pendence on the wave vector, which might be useful to 
characterize dipolar interactions in experimental realiza- 
tions of spin chains. Using the numerical density-matrix 
renormalization group method, we have presented quan- 
titatively accurate results for the magnetization curve 
and the location of the critical point for the case of 
spin S = 1/2, finding a strong reduction in the value 
of h c in comparison to the classical limit (S — > oo). Fi- 
nally, analyzing an effective low-energy field theory for 
our model by means of the functional renormalization 
group method, we have pointed out the emergence of 
a scale-invariant excitation spectrum in the vicinity of 
the critical point, where the logarithmic correction to 
the magnon velocity becomes an irrelevant perturbation, 
thus confirming the consistency of our numerical DMRG 
results with the expected universality class. 



ACKNOWLEDGMENTS 

This work was financially supported by the DFG via 
SFB/TRR49 and FOR 723. 



6 



1 M. A. Baranov, Phys. Rept. 464, 71 (2008). 

2 J. Stuhler, A. Griesmaier, T. Koch, M. Fattori, T. Pfau, S. 
Giovanazzi, P. Pedri, and L. Santos, Phys. Rev. Lett. 95, 
150406 (2005). 

3 T. Lahaye, T. Koch, B. Frohlich, M. Fattori, J. Metz, A. 
Griesmaier, S. Giovanazzi, and T. Pfau, Nature 448, 672 
(2007). 

4 M. Lu, N. Q. Burdick, S. H. Youn, and B. L. Lev, 
Phys. Rev. Lett. 107, 190401 (2011). 

5 K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, B. 
Neycnhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, D. 
S. Jin, and J. Ye, Science 322, 231 (2008). 

6 A. Friedenauer, H. Schmitz, J. Gliickert, D. Porras, and T. 
Schatz, Nature Phys. 4, 757 (2008). 

7 K. Kim, M.-S. Chang, S. Korenblit, R. Islam, E. E. Ed- 
wards, J. K. Freericks, G.-D. Lin, L.-M. Duan, and C. 
Monroe, Nature 465, 590 (2010). 

8 A. I. Akhiezer, J. Phys. USSR, 10, 217, (1946). 

9 S. V. Maleev, Zh. Eksp. Teor. Fiz. 70, 2374 (1976) [Sov. 
Phys. JETP 43, 1240 (1976)]. 

10 A. V. Syromyatnikov, Phys. Rev. B 74, 014435 (2006); 
ibid. 77, 144433 (2008). 

11 J. P. Kestner, B. Wang, J. D. Sau, and S. Das Sarma, 
Phys. Rev. B 83, 174409 (2011). 

12 A. V. Gorshkov, S. R. Manmana, G. Chen, J. Ye, E. Dern- 
ier, M. D. Lukin, and A. M. Rey, Phys. Rev. Lett. 107, 



115301 (2011); A. V. Gorshkov, S. R. Manmana, G. Chen, 
E. Dernier, M. D. Lukin, and A. M. Rey, Phys. Rev. A 84, 
033619 (2011). 

E. G. Dalla Torre, E. Berg, and E. Altman, 
Phys. Rev. Lett. 97, 260401 (2006). 

H. S. Singh, B. B. Singh, R. K. Misra, G. P. Gupta, and 
K. C. Lai, J. Phys. C: Solid State Phys. 19, 4749 (1986). 
L. Lewin, Polylogamthms and Associated Functions, 
(North-Holland, New York, 1981). 

Q. P. Li and S. Das Sarma, Phys. Rev. B 43, 11768 (1991). 

H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 

H. Inoue and K. Nomura, J. Phys. A: Math. Gen. 39, 2161 

(2006). 

S. R. White, Phys. Rev. Lett. 69, 2863 (1992); 
Phys. Rev. B 48, 10345 (1993); R. M. Noack and S. R. 
White, in Density-Matrix Renormalization, edited by I. 
Peschel, X. Wang, M. Kaulke, and K. Hallberg, (Springer, 
Berlin, 1999); U. Schollwock, Rev. Mod. Phys. 77, 259 

(2005) . 

For a recent review see, for example, P. Kopietz, L. Bar- 
tosch, and F. Schiitz, Introduction to the Functional Renor- 
malization Group, (Springer, Berlin, 2010). 
N. Hasselmann and P. Kopietz, Europhys. Lett. 74, 1067 

(2006) . 

D. F. Litim, Phys. Rev. D 64, 105007 (2001). 
M. E. Fisher, Rev. Mod. Phys. 46, 597 (1974); K. G. Wil- 
son, Rev. Mod. Phys. 47, 773 (1975). 



